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'-^ ; Abstract 

Z^ \ New definitions of Q-conditional symmetry for systems of PDEs are presented, whicii 

generalize the standard notation of non-classical (conditional) symmetry. It is shown 

' (— i '' that different types of Q-conditional symmetry of a system generate a hierarchy of condi- 

Q^, tional symmetry operators. A class of two-component nonlinear RD systems is examined 

(— I \ to demonstrate applicability of the definitions proposed and it is shown when different 

definitions of Q-conditional symmetry lead to the same operators. 

PACS numbers: 02.20.-a, 02.30.Jr, 05.45.-a 

1. Introduction 



C^ 



(N 

^ \ Since 1952, when A.C Turing pubhshed the remarkable paper [1], in which a revolutionary 

^r^ \ idea about mechanism of morphogenesis (the development of structures in an organism during 

t^^ ' its life) has been proposed, nonlinear reaction-diffusion systems have been extensively studied 



>< 



by means of different mathematical methods, including group-theoretical methods (see [21 |3l 
m |5] and the papers cited therein) . The main attention was paid to investigation of the two- 



o, 

O . component RD systems of the form 



Ut = [D\U)U^], + F{U,V), 
Vt = [D\V)V,], + G{U,V) 



(1) 



c^ ■ where U = U{t, x) and V = V{t, x) are two unknown functions representing the densities of 

populations (cells, chemicals etc.), F{U,V) and G{U,V) are two given functions describing 
interaction between them and environment, the functions D^{U) and D'^iV) are the relevant 
diffusivities (hereafter they are assumed positive functions) and the subscripts t and x denote 
differentiation with respect to these variables. Note that nonlinear system ([1]) generalizes many 
well-known nonlinear second-order models used to describe various processes in physics [7|, 
biology jSl 13 and ecology [TU] • 

At present, one can claim that all possible Lie symmetries of ([1]) with the constant diffu- 
sivities are completely described in [21 13], while in [6l |5] it has been done for the non-constant 
diffusivities. However, the problem of construction of Q-conditional symmetries (non-classical 
symmetries) for ([T]) is still not solved. Moreover, to our best knowledge, there are only a few 
papers devoted to the search of conditional symmetries of evolution systems, which contain RD 
systems as a subclass [11] -[H]. Notably, some general results about Q-conditional symmetries 
of systems with power diffusivities of the form 

Ut = {U''U,), + F{U,V), , . 

Vt = {V'V,)^ + G{U,V) ^^> 



have been obtained in the recent paper [Tl|. However, the results obtained in [13] cannot 
be adopted for any system of the form ([T]). Moreover, the authors clearly indicated that two 
interesting cases, I = k = and I = k = —1, wasn't studied therein. 

It should be noted that, starting from the pioneering work [15], many papers were devoted to 
the construction of such symmetries for the scalar non-linear reaction-diffusion (RD) equations 
of the form [IS]-[2Z] (the reader may find more references in the recent book [2S] ) 

Ut = [DiU)U,l + F{U). (3) 

and ([3]) with the convective term B{U)Ux (here B{U),D{U) and F{U) are arbitrary smooth 
functions). It is well-known that conditional symmetries can be applied for finding exact 
solutions of the relevant equations, which are not obtainable by the classical Lie method. 
Moreover the solutions obtained in such a way may have a physical or biological interpretation 
(see, e.g., examples in [221 123 ESI EZ])- Thus, the same should be expected in the case of RD 
systems. 

In the present paper, we shall consider a general multi-component system of evolution 
PDEs and its particular case, system ([1]). The paper is organized as follows. In the second 
section, we present new definitions of Q-conditional symmetry, which create a hierarchy of 
conditional symmetry operators. In the third section, the definitions are used to construct 
the systems of determining equations(DEs) for the RD system ([1]). An analysis of the derived 
systems is carried out, particularly, an example is presented, which illustrates different types 
of Q-conditional symmetries. Finally, we present some conclusions in the last section. 

2. Definitions of conditional symmetry for systems of PDEs 

Here we present new definitions of Q-conditional symmetry which naturally arise for systems 
of PDEs. To avoid possible difficulties that can occur in the case of arbitrary system of PDEs, 
we restrict ourself on systems of evolution PDEs (the RD system ([1]) is a typical example). 

Consider a system of m evolution equations {m > 2) with 2 independent {t, x) and m 
dependent u = {ui, U2, ■ ■ ■ , tim) variables. Let us assume that the A;-th order (A; > 2) equations 
of evolution type 

ui = F' {t, X, u, Ux, ..., 4^')) , z = 1, 2, . . . , m, (4) 

are defined on a domain Vl C R^ of independent variables t and x. Hereafter F* are smooth 
functions of the corresponding variables, the subscripts t and x denote differentiation with 

respect to these variables, u^J^ = ^= (|^, • • • , ^) , j = 1, 2, . . . , A;,. 

It is well-known that to find Lie invariance operators, one needs to consider system (jl]) as 
the manifold Ai = {Si = 0, S2 = 0, . . . , S^ = 0}, where 

S, = ui- F' [t, X, u, Ux, . . . , 4'=')) = 0, i = 1, 2, . . . , m, (5) 

in the prolonged space of the variables: 

t,x,u,u,. 



Here k = maxjfcj, i = 1, . . . , m} and the symbol u (^j = 1^2, ... ,k) denotes totalities of the 
j-th order derivatives w.r.t. the variables t and x. 

According to the definition, system (^ is invariant (in Lie sense !) under the transformations 
generated by the infinitesimal operator 

Q = e{t,x,u)dt + e{t,x,u)d^ + r]\t,x,u)du, + . . . + ii"'{t,x,u)d^^, (6) 

if the following invariance conditions are satisfied: 

QSi = Q{ui-F'{t,x,u,u,,...,ui'^^)) =0, 2 = 1,2,. ..,m. (7) 

k k M 

The operator ^ is the k-th order prolongation of the operator Q and its coefficients are expressed 
via the functions ^*^, ^^, r^^, . . . , r]*" by well-known formulae (see, e.g., |29]-|30]). 

The crucial idea used for introducing the notion of Q-conditional symmetry (non-classical 
symmetry) is to extend the manifold A^, namely: the operator Q is used to create a such 
extension. It can be noted that there are several different possibilities to realize this idea in the 
case of systems containing m PDEs. In fact, different definitions can be formulated for such 
systems depending on the number of complementary conditions generated by the operator Q. 
Definition 1. Operator ([6]) is called the Q-conditional symmetry of the first type for an 
evolution system of the form (jl]) if the following invariance conditions are satisfied: 

QS-, = Q(iij-F*(t,a;, «,«,,..., 4'=^))) =0, z = l,2,...,m, (8) 

where the manifold A^i is one {Si = 0, 5*2 = 0, . . . , Sm = 0, Q{ui^) =0} with a fixed number 

H (1 < ^1 < ''^)- 

Definition 2. Operator (E]) is called the Q-conditional symmetry of the p-th type for an 
evolution system of the form (|4]) if the following invariance conditions are satisfied: 

QSi = Q{ul- F' {t,x,u,u,, . . . ,4'=^))) =0, 1 = 1,2,. ..,m, (9) 

where the manifold A^p is one {Si = 0, 52 = 0, . . . , S"™ = 0, Q{ui^) = 0, . . . , Q{ui^) = 0} with 
any given numbers ii, . . . ,ip{l < p < ip < m). 
Definition 3. Operator ([6]) is called the Q-conditional symmetry (non-classical symmetry) 
for an evolution system of the form (jl]) if the following invariance conditions are satisfied: 

QSi = QU^^-F'{t,x,u,u,,...,u'i^'^)) =0, z = l,2,...,m, (10) 

" "^ Mm 

where the manifold M-m is one {Si = 0, 5*2 = 0, . . . , Sm = 0, Q{ui) = 0, . . . , Q{u.m) = 0}. 

Obviously, all three definitions coincide in the case of m = 1, i.e. a single evolution equation. 
If m > 1 then one obtains a hierarchy of conditional symmetry operators. It is easily seen that 
Al C A^i C A^p C M-m^ hence, each Lie symmetry is automatically a Q- conditional symmetry 
of the first and p-ih types, while Q-conditional symmetry of the first type is that of the p-th 
type. From the formal point of view is enough to find all the Q-conditional symmetry (non- 
classical symmetry) operators. Having the full list of Q-conditional symmetries one may simply 
check which of them is Lie symmetry or/and Q-conditional symmetry of the p-th type. 



On the other hand, to construct any Q-conditional symmetry for a system of PDEs, one 
needs to solve another nonhnear system, so called system of DEs, which usually is much more 
complicated and cumbersome. This problem arises even in the case of single linear PDE and it 
was the reason why G.Bluman and J.Cole in their well-known work |15] were unable to describe 
all Q-conditional symmetries in explicit form even for the linear heat equation. Thus, all three 
definitions are important from theoretical and practical point of view. Notably, definition 1 
should be more applicable for solving systems of DEs when one examines a multi-component 
systems containing three and more PDEs. 

It should be stressed that Definition 3 was only used in papers [II]-[ll] devoted to the search 
Q- conditional symmetries for the systems of PDEs. Moreover, to our best knowledge, nobody 
has noted that several definitions producing a hierarchy of conditional symmetry operators can 
be defined for systems of PDEs. 

3. Conditional symmetries for the RD systems 

Consider the RD system ([1]). According to the definitions presented above, two types 
of conditional symmetry operators can be derived: Q-conditional symmetry of the first type 
and Q-conditional symmetry of the second type. The second type coincides with standard 
non-classical symmetry. It turns out that systems of DEs corresponding to both types have 
essentially different structure. 

First of all, the RD system ([T]) can be simplified if one applies the Kirchhoff substitution 



u 



D\U)dU, v= D\V)dV, (11^ 



where u{t, x) and v{t, x) are new unknown functions. Hereafter we assume that there exist 
inverse functions to those arising in right-hand-sides of ( ITTj) . Substituting (ITT|) into ([T]), one 
obtains 

u^x = d^iu)ut + C'^{u,v), ,^2) 

Vxx = d'^{v)vt + C^{u,v), 

where the functions d^, d'^ and C^, C^ are uniquely defined via D^, D^ and F, G, respectively. 
Consider a conditional symmetry operator of system ([T]) 

Q, = C°(t, X, U, V)dt + C'(t, X, U, V)d, + vlit, X, U, V)du + rjlit, x, U, V)dv, (13) 

where ^^,^1,1^1 and rj"^ are to-be-found functions. The Kirchhoff substitution transforms f lT^ 
to the same form 

Q = ^°(t, X, u, v)dt + ^^(t, X, u, v)d,j; + r]^{t, x, u, v)du + rfit, x, u, v)dy. (14) 

where the operator coefficients are uniquely expressed via those with stars using formulae (TTT]) . 
Let us use definition 1 to construct the system of DEs for finding Q-conditional symmetry 
operators of the first type. According to the definition the following invariance conditions must 
be satisfied: 



Q{d\u)ut + C\u,v)-u,,) 
Q{d\v)vt + C\u,v)-v,,) 



= 0, 

^1 (15) 

= 0, 

Ml 



where the manifold 

Ml = {Si = 0,82 = 0, Q{u) = ^\t, X, u, v)ut + ^\t, X, u, v)u^ - r]\t, x, u, v) = 0}. (16) 

Note that the condition Q{v) = instead of Q{u) = can be also used, however, it will lead to 
such conditional symmetry operators, which are obtainable from those generated by the invari- 
ance conditions flT5l) - flT6l) . In fact, the simple renaming u ^ v,v ^ u (trivial discrete transfor- 
mations) and the corresponding coefficients' renaming preserves the form flT2l) and transform 
the invariance criteria flTSj) into that with the condition Q{v) = because arbitrariness of the 
functions d'' and C^, A; = 1, 2. 

Remark 1. In the case of system flT^ with the fixed functions d'^ and C'^, k = 1,2 the 
discrete invariance u ^ v,v ^ u can be broken, so that definition 1 in two cases should be 
examined. 

Now we apply the rather standard procedure for obtaining system of DEs, using the in- 
variance conditions (1151) . From the formal point of view, the procedure is the same as for 
Lie symmetry search, however, three (not two !) different derivatives, say Uxx,Vxx,Ut, can be 
excluded using the manifold A^i. After rather cumbersome calculations, one arrives at the 
nonlinear system of DEs 

-*-/ Sx ><u ^v S-u Sti U, 

^) Vv ~ Vuu ~ Vuu ~ Vvv ~ Vuv ~ U) 

3)evlid'-d') + 2evL = 0, 

4) iae - eel - 2eQd' - iWdi - 2^%^ + eei^ = o, 

5) (2^^ - e)d' + V'dl = 0, 
6)C}d' + 2rjl-CL = 0, 

7) $v'di + ivl + 2a$ - e$)d' - VL + v'Ci + v'Cl + (2^ - vl)C' = 0, 

8) ivf + 1.vl)d' - '^Vld' - vL + V'CI + r^'Cl - r^C' + (2^^ - ^l)C^ = 0, 



(17) 



if ^° 7^ and d^{u) ^ d'^{v) (the special cases ^° = and d^{u) ^ d'^{v) = const can be treated 
in the similar way, however, the results omitted here to avoid cumbersome formulae). 



An analysis of system (11) shows that finding of Q- conditional symmetry of the first type 
cannot be reduced to the operators ( [7^ with ^ = ^- This is in contradiction to the well-known 
fact occurring in the case of single evolution equations. Indeed, if one assumes that system flTTl) 
is locally equivalent to that with ^° = 1 then equation 5) takes the form 

2ej' + v'dl = 0, 

therefore the restriction ,^^ = is obtained provided c?^ = const ^ 0. On the other hand, there 
are several RD systems, which are invariant under the Lie symmetry operators of the form (see 
[2], table 3) 

2tdt + xdx + aiudu + a2vdy, 



(here ai and a2 are constants), i.e., those with ^^ 7^ 0. Obviously, such operators will be lost 
if one assumes ^° = 1 in the very beginning. 

Now we apply definition 3 to construct the system of DEs for finding Q-conditional sym- 
metry (non-classical symmetry) operators. Thus, the invariance conditions 



Q{d\v)vt + C\u,v)-v,,) 



M2 
M2 



0, 







(18) 



must be satisfied. Here the manifold M.2 = {'S'l = 0, 5*2 = 0, Q{u) = 0, Q{v) = 0}, so that four 
different derivatives, say u^x, Vxx, Ut, Vt, can be excluded in this case. Finally, the system of DEs 



/ ^X ?U ^D ><UU ^VU ^VV 

2) vl = vL = 0, 

3) 2|ie^di + vl^ - 2^1 = 0, 

5) |^C'(t/i + rf2) + 2771,-2^^ = 0, 
6)|ie(rfi + rf2) + 2r/i-2eL = 0, 



0, 



7) 1^771(^1 - <f) - 2'^e.d' + 2r,l - 2CC' = 0, 

8) $vlid' - d') - 2^^eud' + 2r/l - 2^0' = 0, 

9) {$^f + 2|^e - ^^ - $a - 4ex)d' - $v'dl+ 

+^^ad' - Hu + exx + ^iic' + ac' = o, 

10) (f^e? + 2$a - ^i - ^a - 2$ex)d' - $v'di+ 

+^^eud' - 2vi + exx + 3CC2 + ac = o, 

11) ivl - $^? + $VI + 2$ex)d' + $v'di - $vld' - tx 

WCi + ri^Cl + (2ei. - r^DC - vlC = 0, 

12) iv^ - $et + $vl + 2$ex)d' + 'i^v'dl - ^^vld' - vlx 
WCl + r^'Cl - vlC' + {2^ - V'.)C' = 



(19) 



is obtained if ^° 7^ (the special case ^° = can be treated in the quite similar way). 

An analysis of system / fTPj) shows that finding Q-conditional symmetry operators for the 
RD system (Cp can be reduced to the operators (fl^ with ^° = 1 provided ^ ¥" ^- 1^ ia-ct, the 
substitution 



e* 



£0' '* 



k = 1,2 



(20) 



reduces system (IT^ to the form (the stars next to the functions ^ and r]'' are skipped): 

-Lj <,UU <,VV SMf ") 

2) ??™ = vL = 0, 

3) 2ee«di + r]i^ - 2U = 0, 

4) 2^^,(P + 7/2^ - 2U = 0, 
5)ee.(t^' + rf') + 2?7i-2e.. = 0, 
Q)iUd^ + (P) + 2ril,-2U = Q, 

7) e?7.Hrf' - rf') - '^ivV'd^ + 2r7L - '^i.C^ = 0, 

8) e^^(^2 _ ^1) _ 2^^^2^2 + 2r^2^ - 2e„C2 = 0, (21) 

9) (2^^771 - 6 - e.^' - 2eex)rf' - e^'<+ 

+e.r72(/2 - 2r^l^ + e... + 3e„Ci + ^.C^ = 0, 

10) (2^,7/2 - e* - e^r^i - 2ii,)<P - ir,^dl+ 

HuTl^d^ - 27]l + ^xx + 3e.C2 + e«Ci = 0, 

11) ivl + rj'vl + 2i,r^^)d^ + {r^'fdl - vWj' ' vL 

WCl + v'Cl + (2e. - ril)C^ - vlC' = 0, 

12) iv? + vWu + 2^.v')d' + iv^dl - vWud' - vL 

Wcl + v'cl - vlc' + (2e. - V'.)C' = 

and ^° is an arbitrary smooth function of the time variable. Now we reahze that system (12T1) 
is nothing else but the system of DEs to find the Q-conditional symmetry operators (TT4|) with 
^° = 1. Thus, to find all Q-conditional symmetries of the RD system ([T]) one needs to solve 
only the particular case of system (1191) . i.e., system (I2ip . 

Remark 2. System fl2Tl) with d^{u) = u^,d'^{v) = w" coincides with the system of DEs 
obtained and analyzed in the recent paper [H]. 

Of course, system of DEs flT7|) can be derived from system flT9l) because each Q-conditional 
symmetry of the first type is automatically a Q-conditional symmetry (but not vice versa !). 
However, system fITIJl) is much complicated than fITTI) and its solving is a difficult problem even 
in particular cases, e.g., d^{u) = u"^,d'^{v) = v"' [H]. Thus, we believe that Q-conditional 
symmetries of the first type can be found much easier for many nonlinear RD systems ([1]) with 
correctly-specified coefficients, arising in applications. 

Now we present an example, which highlights when Q-conditional symmetry (non-classical 
symmetry) is or not a Q-conditional symmetry of the ffist type. 
Example. Consider system (I19p assuming 

e = Vl = vl = 0- (22) 

In this case, the system can be immediately simplified and one obtains only two equations 

ivl - $C?)d' + $v'di - vL + ri'Cl + v'Cl - r^lC' = 0, 

{ri't - ''iet)d' + $^r,^dl - r^l + r^'C^ + r,^Cl - ^C^ = ^ ^ 



to find the functions 

r]^ = r\t)u + p^{t,x), (24) 

7^2 = r'^(t)v +p'^{t,x), 

wliere tlie functions in rigfit-liand-sides should be determined. 
System (fT7|) with restrictions (!22|) takes the form 

ivl - $C?)d' + $v'di - vL + V'CI + ri^Cl - viC' = 0, 

v^d' - vL + v'cl + v'c'. - nlc^ = O, (25) 

Cy = V'dl 



where the unknown functions C,^ and r]^,k = 1,2 must be of the form (^ 

Now one notes that systems fl23l) and fl25l) are equivalent if rj'^dl = (we remind the 
reader that substitution fl20|) reduces system flTOl) to that with ,^° = 1). Thus, we arrive 
at the statement: each Q-conditional symmetry operator [14\ ) of the RD system [W^) with 
(P{v) = const (or d^{u) = const) is equivalent (up to multiplier c{t) ^ 0) to a Q-conditional 
symmetry operator of the first type if restrictions [2^) take place. Moreover, the statement is 
valid for all systems of the form [W\) under the additional restriction rj^ = (or rj^ = 0). 

Taking into account the Kirchhoff substitution, the same statement can be formulated for 
RD system ([T]) because substitution ( ITTl) preserves restrictions (1221) . 

Finally, we analyze briefly the RD system ([2]) (the cases I = k = and I = k = —1 
are excluded) using paper table 1 [Hj, which presents all possible Q-conditional symmetry 
(non-classical symmetry) operators of ([2]) under restrictions fl22|) . There are five different cases 
according to theorem 1 p^ . Using the statement obtained we conclude that the Q-conditional 
symmetry operator arising in the first case of table 1 is simultaneously a Q-conditional sym- 
metry operator of the first type while the operator arising in the third case is not because 
I = k = — I 7^ and rj^rf ^ 0. Other three operators of Q-conditional symmetry from Table 
1 are those of Q-conditional symmetry of the first type only under additional restrictions on 
coefficients of the corresponding RD systems. These restrictions can be easily derived using 
the statement formulated above. For example, the Q-conditional symmetry operator listed the 
second case is that of the first type under the restrictions either Ai = or / = (if Ai = / = 
then it is the Lie symmetry operator of the corresponding RD system). 



Table 1. Q-conditional symmetries of the RD system ( I12j) 



RD systems of the form ([T]) 



Q-conditional operators 



Restrictions 



Vt = (V-^V^), - 2XVi + .g(C/'=+i) 



dt + 2p{x)V2dv 



Pxx = (p)^ + Ap, 






dt + \iU-''du+ 
+X2V-^dv 



_ Ai(fc+1) ^ / f. 

Af + Z V 0. 



C/t = iU-2U,), - 2AC/2 + /(C/2 _ y 2), 
\4 = (y~^Vx)x - 2XVi + g{ui - V^) 



dt+ 

-2p{x){uidu + vidv) 



Pxx = [pf + Ap, 






5t + AiC/-^%+ 
+X2{V ~ \:iV-^)dv 



cxpa''+^ 



Ai(fc + 1) 7 A2 7^ 0, 

either A^ + A§ 7^ or 
Ai + /c V or A? + ; V 



f/t = (f/*=C/,). + (C/'=+i - Ai)(/(l.) + A2C/-'=) 
Vt = iV'V,), + (V^'+i - \^){9{lo) + \iV-') 



dt + \2{U-\iU~^)du + 

+\i{V - \:iV-^)dv 



U} 



2{k + l) 1 A2A4 7^ 0, 



A2(fc + 1) 

either Af + A§ 7^ 0, 
or A§ + fc V or Af + / V 



4. Conclusions 

In this paper, new definitions of Q-conditional symmetry for systems of PDEs are presented, 
which generalize the standard notation of non-classical (conditional) symmetry. It is shown 
that different types of Q-conditional symmetry generate a hierarchy of conditional symmetry 
operators. 

Systems of DEs to find Q-conditional symmetries of two types for the nonlinear RD system 
(II]) are constructed. The case of conditional invariance under operator (TT3|) with the coefficients 
^^ = 0, rjlitjXjU) and r]l{t,x,V) is analyzed in details. Using the recent paper [H], we 
established that there are the nonlinear RD systems of five types, which possess Q-conditional 
symmetry (non-classical symmetry) operators satisfying definition 3. However, only one of them 
is simultaneously a Q-conditional symmetry operator of the first type according to definition 
1. 

It is well-known that conditional symmetries can be applied for finding exact solutions of 
the relevant equations, which are not obtainable by the classical Lie method. In a forthcoming 
paper, we are going to find and apply conditional symmetry operators of different types for 
constructing new exact solutions of the diffusive Lotka-Volterra system, which is one of the 
most popular RD systems. 

Finally, we note that definitions 1-3 can be straightforwardly extended on an arbitrary m- 
component system of PDEs presented in a 'canonical' form (the system has a simplest form 
and there are no any non-trivial differential consequences). However, new difficulties may arise 
because differential consequences of additional conditions 



Q(mi) =0,...,Q(m„) = 0, 



(26) 



generated by the conditional symmetry operator (jS]). In fact, there are examples of single 
hyperbolic PDEs when the set of Q-conditional symmetry operators can be extended if one takes 
into account such differential consequences. For example, this occurs in the case of nonlinear 
hyperbolic equation uu = uu^x- The reader may easily check that there is the much wider set of 
Q-conditional symmetry operators in [16] (Chapter 5) than it was found in [31] (Supplement 7) 
neglecting differential consequences. It is the reason why there is the definition of Q- conditional 
invariance (non-classical invariance) for a fc-th order PDE, which requires to take into account 
all differential consequences of additional condition up to the order k (see, e.g., [32]). From this 
point of view, definitions 1-3 can be formally extended by using differential consequences of (126!) 
up to the order k = m.ax{ki, i = 1, . . . , m}. However, we believe that such extensions produce 
much more cumbersome formulae but don't lead to any new operators for evolution system (jlj) 
because all the above mentioned differential consequences contain the mixed and/or high-order 
time derivatives, Utx,uu,Uttx, ■ ■ ■, which are absent in system (jl]). For example, if one applies 
the 'amended' definition 1 with Ml = {Si = 0, 5*2 = 0, Q{u) = 0, dtQ{u) = 0, dxQ{u) = 0} 
(instead of flT6|) ) to find Q-conditional symmetry operators of the first type for the RD system 
flT2|) then the system of DEs fITTj) is exactly obtained but nothing more. 
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^ • 1. Introduction 

y~^ I Since 1952, when A.C Turing published the remarkable paper [T], in which a revolutionary 

O I idea about mechanism of morphogenesis (the development of structures in an organism during 

^-^ • its life) has been proposed, nonlinear reaction-diffusion systems have been extensively studied 

by means of different mathematical methods, including group-theoretical methods (see [21 |3l 
m |5] and the papers cited therein) . The main attention was paid to investigation of the two- 
rS I component RD systems of the form 



Ut = [D\U)U,l + F{U,V) 
Vt = [D\V)V,l + G{U,V) 



2.mT.i , r^nr... (1) 



where U = U{t, x) and V = V{t, x) are two unknown functions representing the densities of 
populations (cells, chemicals etc.), F{U,V) and G{U,V) are two given functions describing 
interaction between them and environment, the functions D^(U) and D'^iV) are the relevant 
diffusivities (hereafter they are assumed positive functions) and the subscripts t and x denote 
differentiation with respect to these variables. Note that nonlinear system ([T]) generalizes many 
well-known nonlinear second-order models used to describe various processes in physics ^, 
biology [TJ [8] and ecology [9] . 

At present, one can claim that all possible Lie symmetries of ([1]) with the constant diffusiv- 
ities are completely described in [2l[3], while it has been done for the non-constant diffusivities 
in [S]. However, the problem of construction of Q-conditional symmetries (non-classical sym- 
metries) for ([T]) is still not solved. Moreover, to our best knowledge, there are only a few papers 



devoted to the search of conditional symmetries of evolution systems, which contain RD sys- 
tems as a subclass [TU] ^[13] • Notably, some general results about Q-conditional symmetries of 
systems with power diffusivities of the form 






(2) 



have been obtained in the recent paper |T3]. However, the results obtained in [13] cannot 
be adopted for any system of the form ([1]). Moreover, the authors clearly indicated that two 
interesting cases, / = A; = and I = k = —1, wasn't studied therein. 

It should be noted that, starting from the pioneering work [13], many papers were devoted to 
the construction of such symmetries for the sca/ar nonlinear RD equations of the form [T5] |23] 
(the reader may find more references in the recent book |24j ) 

Ut = [DiU)U.l + FiU). (3) 

and ([3]) with the convective term B{U)Ux (here B{U),D{U) and F{U) are arbitrary smooth 
functions). It is well-known that conditional symmetries can be applied for finding exact 
solutions of the relevant equations, which are not obtainable by the classical Lie method. 
Moreover the solutions obtained in such a way may have a physical or biological interpretation 
(see, e.g., examples in [191 ISSl [23]). Thus, the same should be expected in the case of RD 
systems. 

In the present paper, we shall consider a general multi-component system of evolution PDEs 
and its particular case, system ([T]). The paper is organized as follows. In the second section, 
we present new definitions of Q- conditional symmetry, which create a hierarchy of conditional 
symmetry operators. In the third section, the definitions are used to construct the systems 
of determining equations(DEs) for the RD system ([1]). An analysis of the derived systems 
is carried out, particularly, an example is presented, which illustrates different types of Q- 
conditional symmetries. In the fourth section, the the definition of Q-conditional symmetry is 
applied to find new exact solutions of the classical Lotka-Volterra system with diffusive terms. 
Finally, we present some conclusions in the last section. 

2. Definitions of conditional symmetry for systems of PDEs 

Here we present new definitions of Q-conditional symmetry, which naturally arise for systems 
of PDEs. To avoid possible difficulties that can occur in the case of arbitrary system of PDEs, 
we restrict ourself on systems of evolution PDEs (the RD system ([1]) is a typical example). 

Consider a system of m evolution equations (m > 2) with 2 independent (t, x) and m 
dependent u = [ui, U2, ■ ■ ■ , u^) variables. Let us assume that the A;-th order [k > 2) equations 
of evolution type 

4 = F' {t, X, u, Ux, ..., 4^')) , z = 1, 2, . . . , m, (4) 

are defined on a domain n C R^ of independent variables t and x. Hereafter F* are smooth 
functions of the corresponding variables, the subscripts t and x denote differentiation with 
respect to these variables, wj = ff , and ui'^ = B= (&' • • • ' ^) , i = 1, 2, • • • , A;,. 



It is well-known that to find Lie invariance operators, one needs to consider system (jl]) as 
the manifold A^ = {Si = 0,82 = 0, ..., Sm = 0}, where 

Si = ul- F' (t, X, -u, Ux, . . . , u^^'^) = 0, i = 1, 2, . . . , m, (5) 

in the prolonged space of the variables: 

t X u u u 

Here k = maxjfcj, i = 1, . . . , m} and the symbol ^ {j = 1,2, ... ,k) denotes totalities of the 
j-th order derivatives w.r.t. the variables t and x. 

According to the definition, system @ is invariant (in Lie sense !) under the transformations 
generated by the infinitesimal operator 

Q = i^{t,x,u)dt + i^{t,x,u)d^ + ri^{t,x,u)du^ + . . . + ?7™(t,x, it)9„„, (6) 

if the following invariance conditions are satisfied: 

QSi = QU-r(t,x,u,u^,...,u^t^)) =Q, 1 = 1,2,. ..,m. (7) 

k k Ai 

The operator Q is the k-th order prolongation of the operator Q and its coefficients are expressed 

k ^^^ ^^^ 

via the functions ^°, ^^, r;^, . . . , r^*" by well-known formulae (see, e.g., [2S]-[2S])- 

The crucial idea used for introducing the notion of Q-conditional symmetry (non-classical 
symmetry) is to change the manifold Ai, namely: the operator Q is used to reduce A^. It is 
important to note that there are several different possibilities to realize this idea in the case of 
systems containing m PDEs. In fact, different definitions can be formulated for such systems 
depending on the number of complementary conditions generated by the operator Q. 

Definition 1. Operator (|6]) is called the Q-conditional symmetry of the first type for an 
evolution system of the form (jl]) if the following invariance conditions are satisfied: 

Qs, = Q{ul-F'{t,x,u,u,,...,ui'^^)) =0, 1 = 1,2,. ..,m, (8) 

where the manifold A^i is one {Si = 0, 52 = 0, . . . , S*™ = 0, Q{ui^) = 0} with a fixed number 
ii {1 < ii < m). 

Definition 2. Operator (E]) is called the Q-conditional symmetry of the p-th type for an 
evolution system of the form (^ if the following invariance conditions are satisfied: 

QSi = Q{ul-F'{t,x,u,u„...,ui'^^)) =0, 1 = 1,2,. ..,m, (9) 

where the manifold A^p is one {Si = 0, S'2 = 0, . . . , S'm = 0, Q{uiJ = 0, . . . , Q{ui^) = 0} with 
any given numbers ii, ■ ■ ■ ,ip{l < p < ip <m). 

Definition 3. Operator ([6]) is called the Q-conditional symmetry (non-classical symmetry) 
for an evolution system of the form (|1]) if the following invariance conditions are satisfied: 

QS, = Q {e^~ F' {t,x,u,u,, . . . ,u^^^'>)) =0, 1 = 1,2,. ..,m, (10) 

" '^ Mm 



where the manifold A^m is one {Si = 0, S2 = 0, . . . , Sm = 0, Q{ui) = 0, . . . , Q{um) = 0}. 

Obviously, all three definitions coincide in the case oim = 1, i.e. a single evolution equation. 
If m > 1 then one obtains a hierarchy of conditional symmetry operators. It is easily seen that 
Aim C. J^p C M.I C A^, hence, each Lie symmetry is automatically a Q-conditional symmetry 
of the first and p-th types, while Q-conditional symmetry of the first type is that of the p-th 
type. From the formal point of view is enough to find all the Q-conditional symmetry (non- 
classical symmetry) operators. Having the full list of Q-conditional symmetries one may simply 
check which of them is Lie symmetry or/and Q-conditional symmetry of the p-th type. 

On the other hand, to construct any Q-conditional symmetry for a system of PDEs, one 
needs to solve another nonlinear system, so called system of DEs, which usually is much more 
complicated and cumbersome. This problem arises even in the case of single linear PDE and it 
was the reason why G.Bluman and J.Cole in their well-known work [H] were unable to describe 
all Q- conditional symmetries in explicit form even for the linear heat equation. Thus, all three 
definitions are important from theoretical and practical point of view. Notably, definition 1 
should be more applicable for solving systems of DEs when one examines a multi-component 
systems containing three and more PDEs. 

It should be stressed that definition 3 was only used in papers [IHl^IIS] devoted to the search 
Q-conditional symmetries for the systems of PDEs. Moreover, to our best knowledge, nobody 
has noted that several definitions producing a hierarchy of conditional symmetry operators can 
be defined for systems of PDEs. 

3. Conditional symmetries for the RD systems 

Consider the RD system ([1]). According to the definitions presented above, two types 
of conditional symmetry operators can be derived: Q-conditional symmetry of the first type 
and Q-conditional symmetry of the second type. The second type coincides with standard 
non-classical symmetry. It turns out that systems of DEs corresponding to both types have 
essentially different structure. 

First of all, the RD system ([1]) can be simplified if one applies the Kirchhoff substitution 



u 



D\U)dU, v= D\V)dV, (11) 



where u{t, x) and v{t, x) are new unknown functions. Hereafter we assume that there exist 
inverse functions to those arising in right-hand-sides of ( ITT]) . Substituting (ITTj) into ([T]), one 
obtains 

Uxx = d^{u)ut + C^(w, v), .^2^ 



where the functions d^, d^ and C^, C^ are uniquely defined via D^, D^ and F, G, respectively. 
Consider a conditional symmetry operator of system ([T]) 

Q* = e°(t, X, [/, V)dt + e* (t, x, U, V)d, + r,l{t, X, U, V)du + r,l{t, x, U, V)dv, (13) 

where ^'^,^l,r]l and rjl are to-be-found functions. The Kirchhoff substitution transforms 0131) 
to the same form 

Q = C.'^{t,x,u,v)dt + ^^{t,x,u,v)da: + T]^ {t, x,u,v)du + rf{t,x,u,v)dy. (14) 



where the operator coefficients are uniquely expressed via those with stars using formulae (ITT]) . 
Let us use definition 1 to construct the system of DEs for finding Q-conditional symmetry 
operators of the first type. According to the definition the following invariance conditions must 
be satisfied: 



Q {d^{u)ut + C^(m, v) - u^^) 



Ml 
Ml 



= 0, 

0, 



(15) 



where the manifold 



Ml = {Si 



0,^2 



0, Q{u) = ^°(t, x, M, v)ut + i^{t^ X, u, v)ux — r]'^{t, x, u, v) = 0}. (16) 



Note that the condition Q{v) = instead of Q{u) = can be also used, however, it will lead 
to such conditional symmetry operators, which are obtainable from those generated by the 
invariance conditions (I15j) - 016j) . In fact, the simple renaming u ^ v,v ^ u (trivial discrete 
transformations) and the corresponding coefficients' renaming preserve the form flT2|) and trans- 
form the invariance criteria flTSl) into that with the condition Q{v) = because of arbitrariness 
of the functions d'' and C'', k = 1,2. 

Remark 1. In the case of system flT^ with the fixed functions d^ and C^, k = 1,2 the 
discrete invariance u -^ v,v ^ u can be broken, so that definition 1 in two cases should be 
examined. 

Now we apply the rather standard procedure for obtaining system of DEs, using the in- 
variance conditions (fT5l) . From the formal point of view, the procedure is the same as for 



Lie symmetry search, however, three (not two !) different derivatives, say u^ 



,Ut, can be 



excluded using the manifold A^i. After rather cumbersome calculations, one arrives at the 
nonlinear system of DEs 



i)4° = ?° = e° = e^ = C = o, 

2) vl = vL = vL = Vvv = vL = 0, 

3)evlid'-d') + 2evL = 0, 

4) iae - eel - 2eQd' - iWdi - 2^%^ + ea. = o, 

5) (2^^ - e)d' + rfdl = 0, 
6)C}d' + 27jl-CL = 0, 

7) '^V'di + ivl + 2e4 - e$)d' - vL + v'Ci + v'Cl + (2^ - vl)C' = 0, 

8) ivf + $vl)d' - $vld' - vL + V'CI + ^'Cl - r^lC^ + (2^^ - r^l)C^ = 0, 



(17) 



if <^° 7^ and d^{u) ^ (^{y) (the special cases ^° = and dl^{u) = d'^iv) = const can be treated 
in the similar way, however, the results omitted here to avoid cumbersome formulae). 

An analysis of system (7?) shows that finding of Q-conditional symmetry of the first type 
cannot be reduced to the operators [l4\ ) with ^ = ^- This is in contradiction to the well-known 



fact occurring in the case of single evolution equations. Indeed, if one assumes that system (1171) 
is locally equivalent to that with (^"^ = 1 then equation 5) takes the form 

2ej' + v'di = 0, 

therefore the restriction ^^ = is obtained provided (P = const ^ 0. On the other hand, there 
are several RD systems, which are invariant under the Lie symmetry operators of the form (see 
[2], tables) 

2tdt + xdx + aiudu + a2vdy, 

(here ai and 0:2 are constants), i.e., those with ^l 7^ 0. Obviously, such operators will be lost 
if one assumes ^° = 1 in the very beginning. 

Now we apply definition 3 to construct the system of DEs for finding Q-conditional sym- 
metry (non-classical symmetry) operators. Thus, the invariance conditions 



Q {d^{u)ut + C^(m, v) - u^^) 
QAd\v)vt + C\u,v)-v,,) 



M2 



M7 



= 0, 




must be satisfied. Here the manifold M.2 = {'S'l = 0, ^2 = 0, Q{u) = 0, Q{v) = 0}, so that four 



different derivatives, say u^ 



,Ut,Vt, can be excluded in this case. Finally, the system of DEs 



/ ^x ^u ><v >>uu ^vu ^vv 

2) vL = vL = 0, 

3) 2|c7^„d + r^^u — 2^^^ = 0, 

4) 2^^e.d' + Vl - 2C = 0, 
5)^^Cid' + d') + 2vL-2Cl = 0, 
Q)$,eu{d' + d') + 2r^l-2a^ = 0, 



0, 



zzitui 



7) l^vlid' - d') - 2^^ild' + 2r^l - 2^0^ = 0, 

8) ^^Vlid' - d^) - 2^^ad' + 2r^L - 2^0' = 0, 

■J) l|OSt + ^|0"4u ^ St ^ Jo<^v ~ 2|o4xJ^ ~ |o^ "'U~T~ 

+|^Crf2 _ 2ryi„ + e,^ + sac' + ilC' = 0, 

10) (f^e? + 2f^c^ - a - $,eu - 4ci)d' - i^v'di+ 

+ |ow« — 2r]^j.^ + ^^^ + 3t,^C + ^u^ =0, 

11) ivl - $a + $VI + 2f^e^)rf^ + $v'di - $vld' 

Wc'^ + v'cl + (2^ - vi)c' - vie = 0, 

12) iv? - $C? + $vl + 2$a)d' + $v'dl - $vld' • 
WCl + v'Cl - vlC' + (2ei - v'.)C' = 



(19) 



Ixx 



Ixx 



is obtained if (^° 7^ (the special case ^"^ = can be treated in the quite similar way). 

An analysis of system / f7P|) shows that finding Q-conditional symmetry operators for the 
RD system (J\) can be reduced to the operators [I4\ ) with ^'^ = 1 provided ^^ ^ 0. In fact, the 
substitution 






^. = ^,v' = i^,k = l,2 (20) 



reduces system ( IT9l) to the form (the stars next to the functions ^ and r]^ are skipped): 

^) Smm ^vv <,uv U, 

2) vl = vL = 0, 

3) 2^^^d' + r,l^ - 2U = 0, 

4) 2^.^' + Til - 2U = 0, 
5)ee.(rf' + rf') + 2?7i-2e.. = 0, 

Q)^Ud' + d^) + 2r,l,-2U = ^. 

7) U{d} - d^) - 2i,ri^d^ + 2r,l^ - 2i,C^ = 0, 

8) iril{d' - d') - 2^^r^'d' + 2r/L - ^CuC' = 0, (21) 

9) (2e„r;i - e* - ^.v' - 2CQd^ - Cv'di+ 

Hvv'd^ - Hu + e^x + 3e«ci + e.c2 = o, 

10) (2^,7^2 _^^_ ^^^1 _ 2^^^)rf2 _ ^^2^2 + 

Huri'd^ - 27]l + e.. + 3e.C2 + e„Ci = 0, 

11) ivl + vY. + 2^.v')d' + iv^fdi - vWJ' - vL 

WC^ + V'CI + (2e. - ril)C^ - vlC' = 0, 

12) iv? + vWu + 2^.v')d' + iv^dl - vWud' - vL 

WCl + ri'Cl - vlC' + (2^. - v'JC' = 

and ^° is an arbitrary smooth function of the time variable. Now we realize that system (12T1) 
is nothing else but the system of DEs to find the Q-conditional symmetry operators (TT^ with 
^° = 1. Thus, to find all Q-conditional symmetries of the RD system ([1]) one needs to solve 
only the particular case of system (1191) . i.e., system (I2ip . 

Remark 2. System fl2Tl) with d^{u) = u"^,d'^{v) = f" coincides with the system of DEs 
obtained and analyzed in the recent paper 



Remark 3. System of DEs to search Lie symmetries of RD systems with non-constant 
diffusivities (see equations (9)-(13) in [5]) differs essentially from systems (fTTl) and (TT9l) . In 
contrary to systems of DEs for Q-conditional symmetries, one for Lie symmetries contains the 
subsystem of equations 

^0 = ^0 = ^1 = pi = r?i = r?^ = 

Su Si5 Sii Stf 'Iv lu 



what essentially simplifies its solving. 

Of course, system of DEs flT7|) can be derived from system flT^ because each Q-conditional 
symmetry of the first type is automatically a Q-conditional symmetry (but not vice versa !). 
However, system ( !T9|) is much complicated than (TT7|) and its solving is a difficult problem even 
in particular cases, e.g., d^{u) = u"^,cP{v) = v"" p^. Thus, we believe that Q-conditional 
symmetries of the first type can be found much easier for many nonlinear RD systems ([1]) with 
correctly-specified coefficients, arising in applications. 

Now we present an example, which highlights when Q-conditional symmetry (non-classical 
symmetry) is or is not of the first type. 

Example 1. Consider system flT9|) assuming 



e' = Vl = vl = 0. (22) 

In this case, the system can be immediately simplified and one obtains only two equations 

ivl - $^>' + $v'di - vL + v'cl + v'cl - vic' = 0, 

to find the functions 

r]^ = r\t)u + p^{t,x), (24) 

7^2 = r'^(t)v +p'^{t,x), 

where the functions in right-hand-sides should be determined. 
System ( ITTl) with restrictions ( l22l) takes the form 

ivl - $C?)d' + '^v'di - vL + V'CI + ri^Cl - viC' = 0, 

V^d' - vL + V'CI + 7^'Cl - nlC^ = 0, (25) 

ad' = V'dl 



where the unknown functions ^° and r]'^,k = 1,2 must be of the form 

Now one notes that systems (1251) and fl2Sl) are equivalent if rj'^dl = (we remind the 
reader that substitution (120|) reduces system (TT9l) to that with ^° = 1). Thus, we arrive 
at the statement: each Q-conditional symmetry operator [T^ of the RD system [W\) with 
d'^{v) = const (or d^{u) = const) is equivalent (up to multiplier c{t) ^ 0) to a Q-conditional 
symmetry operator of the first type if restrictions [2E) take place. Moreover, the statement is 
valid for all systems of the form fT^) under the additional restriction 77^ = (or r/-*^ = Oj. 

Taking into account the Kirchhoff substitution, the same statement can be formulated for 
RD system ([1]) because substitution fITT]) preserves restrictions fl22l) . 

Finally, we analyze briefiy the RD system ([2]) (the cases I = k = and I = k = —1 
are excluded) using paper table 1 [T3|, which presents all possible Q-conditional symmetry 
(non-classical symmetry) operators of ([2]) under restrictions fl22]) . There are five different cases 
according to theorem 1 [13]. Using the statement obtained, we conclude that the Q-conditional 
symmetry operator arising in the first case of table 1 is simultaneously a Q-conditional sym- 
metry operator of the first type while the operator arising in the third case is not because 



/ = fc = — I 7^ and 77^77^ 7^ 0. Other three operators of Q-conditional symmetry from table 
1 are those of Q-conditional symmetry of the first type only under additional restrictions on 
coefficients of the corresponding RD systems. These restrictions can be easily derived using 
the statement formulated above. For example, the Q-conditional symmetry operator listed the 
second case of table 1 is that of the first type under the restrictions either Ai = or / = (if 
Ai = / = then it is the Lie symmetry operator of the corresponding RD system). 

Table 1. Q-conditional symmetries of the RD system ([2]) 



no 



RD systems of the form ([T|) 



Q-conditional operators 



Restrictions 



Vt 



-^{V--2V^),-2\Vi+giU''+^) 



dt + 2pix)V2dv 



Pxx ^ {pY + Ap, 






dt + XiU-^du 
+X2V-^dv 



\\+P^ 0. 



Ut = {u—2U^)^ - 2XU2 + f{u-2 - y 2), 

Vt = {V-^Vx)x - 2XVi + gjui - V^) 



dt+ 
+2p{x)iU^du + vidv) 



Pxx = (p)'^ + Xp, 
P^O 



Ut = {U>'U^).:, + XiU-'' + f{oj), 
Vt - iV'Vx)x + iV'+' - X3)ig{u;) + X2V-') 



dt + XiU-''du+ 

-X2{V - X3V-^)dv 



cxp U 



fc+1 



A2y^0, 



Ai(fc + 1) 7 

(y'+i-As) A2('+i) 
either A^ + A§ ^ or 
Aj + fc V orXl + P^0 



Ut = [U^U^^y, + (C/^+i - Ai)(/(c^) + AaC/-'^) 
Vt - iV'V,), + {V'+^ ~ X3){g{Lo) + A4I/-O 



dt + X2{u-XiU-'')du- 

+Xi{V - X3V-^)dv 



C/^ + ^-Ai 



U+i) 1 A2A4 7^ 0, 

eit/ier Af + A^ 7^ 0, 
or A2 + A:2 ^ or Af + ^2 ^ 



4. Conditional symmetries of the first type for the diffusive Lotka-Volterra system 



Here we consider the diffusive Lotka-Volterra (DLV) system 



Ai^i 

A2ft 






. + u{ai + hiu + Cif), 



(26) 



which is the most common particular case of reaction-diffusion (RD) system (fT2|) . The param- 
eters Xk,ak,bk, and c^ (A; = 1, 2) are assumed to be arbitrary non-zero constants. 

System ([1]) is the simplest generalization of the classical Lotka-Volterra system that takes 
into account the diffusion process for interacting species (see terms u^x and v^x)- Nevertheless 
the classical Lotka-Volterra system was introduced by A.J.Lotka and V.Volterra more than 
80 years ago, its natural generalization fl26|) is still studied because this is one of the most 
important mathematical models. Lie symmetries of fl26|) have been completely described in [29] 
(note those can be extracted from more general results presented in [21 [3]). 

Theorem 1 In the case Ai 7^ A2, DLV system [2D\) is invariant under Q-conditional operators 
of the first type if and only if the corresponding system and Q-conditional symmetries (up to 



local transformations u — )■ hu^ f — )■ cf , hc^ Q ) have the forms 

XiUt = u,j,,j, + u{ai + u + v) , , . 

\2Vt = V:y,:y, + v{a2 + u + v), ai^a2, 

Qi = (Ai - X2)dt + (ai - a2)uidu - d,), (28) 

Q2 = (Ai - X2)dt + (ai - a2)vid, - 9„). (29) 

Each other Q- conditional operator of the first type coincides with the Lie symmetry operator. 

In the case Ai = A2, DLV system j^) is invariant only under such Q-conditional operators 
of the first type, which coincide with the Lie symmetry operators. 

Proof. Proof is based on solving the system of DEs fITTI) with d'^ = A^, C^ = u{ai + biu+civ) 
and C^ = v{a2 + 62M + C2v). After rather simple calculations, one obtains that DLV system 
( 126|) is invariant only under a Q-conditional operator of the first type only under the coefficient 
restrictions 

Ai 7^ A2, 61 = 62 = b, ci = 02 = c. 

In this case the transformations u — )■ bu, f — )■ cf reduce the DLV system to the form (!27|) . The 
corresponding operator takes the form ( !28|) . 

To find operator (!29|) one needs to proceed the same algorithm but with the invariance 
criteria ( !T5|) for the manifold 

Ml = {Si = 0,82 = 0, Qiv) = e°(t, X, u, v)vt + ^\t, X, u, v)v, - ri^it, x, u, v) = 0}, (30) 

i.e., to use definition 1 with ii = 2. 

Remark 4. Theorem 1 gives a complete description of Q-conditional symmetries of the 
first type in explicit form for DLV system fl26|) . It turns out that all possible Q-conditional 
symmetries of the second type cannot be derived in a similar way. In fact, the system of DEs 
f l?I]) with d^ = Afc, C^ = u{ai + biu + civ) and C^ = v{a2 + b2U + C2v) has much more complicated 
structure and can be solved only under some additional assumptions. 

The Q-conditional symmetry operators (12 8 p and (12 9 p can be applied for finding exact solu- 
tions of DLV system fl27|) using the same algorithm as for classical Lie symmetries. Since both 
operators have the same structure we examine only the second one. Thus, the corresponding 
ansatz reads 

u{t,x) = ^i{x)-^2{x)exp{f^^t), 

^;(t,x) = ^2(x)exp(M|t) ^^ ' 

and reduces system (p7j) to the ODE system 

if'l + fl + ai^i = 0, 

// _|_ a2Ai-aiA2 ,„ J_ ,n ,n — n ^ ) 

V2 H A1-A2 ^2 + V'lV'2 - U, 

in 



where (pi{x) and (p2{x) are to be determined functions. 

The general solution of this nonlinear ODE system cannot be found in an explicit form. 
However, we constructed its particular solutions, which lead to interesting solutions of the 
corresponding DLV system. In fact, the first ODE in (1291) possesses the non-zero steady state 
solution Lpi = —tti. Substituting this solution into the second ODE in ( 129|) . one arrives at the 
linear ODE 

^2 - /3Xm = 0, (33) 

where /? = y~^ y^ 0. Thus, using the general solution of fl5^ and ansatz fl^ . we obtain two 
families of exact solutions of the DLV system (^7^: 

u{t,x) = -ai + -^— (Ci exp{y//3Xix) + C2exp{-^/^>^x))e'^\ 
v{t,x) = ^^{Clexp{^/^x) + C2exp{-^/pX^x))e^', 

if ^ > 0, and 

u{t,x) = — ai H ^ (Ci cos{y/—/3Xix) + C2sin(V— /3Aix))e^*, 

v{t,x) = ^(Cicos(v/^:5ATa;) +C2sin(v/^:5Arx))e^*, 

if /? < (hereafter Ci, C2 are arbitrary constants). 

It should be noted that solutions ( !M|) and (135!) cannot be constructed using Lie symmetries. 
In fact, DLV system ( l27l) admits only the trivial Lie symmetry J29], hence, the plane wave 
solutions only can be found by Lie symmetry reductions (the reader may find examples of 
such solutions in [30] and [29]). Obviously, solutions ( 1M1) and ( l35l) possess more complicated 
structures. 

Example 2. Consider solution ( 135|) with Ci = 0. Using the substitution u — )• —bu, v — )■ 
—cv {b > 0, c > 0), one transforms DLV system ( !27l) to the standard system describing the 
competition of two species (see, e.g., [7]): 

XiUt = U:^x + u{ai - bu - cv) , , - 

X2Vt = v^^ + v{a2 -bu- cv) 

and solution (l35ll to the form 

^t, x) = f + j,^^ C2 sm{V^Wix)e^\ 
^('^'^) = (a2~ai)c C2sm{y/-^Xix)e^\ 

where the coefficient restrictions /3 = I^Z'x < 0, ai > 0, 02 > are assumed. 
One notes that this solution satisfies the constant Dirichlet conditions 

x = 0: u = f, v = 0, 

T = . "",,. • ?; = 21 7; = n ^ > 
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in the domain fi = {{t,x) E (0, +00) x (o^-^===\j. Moreover, solution (l37|l has the time 
asymptotic 

(m, v)^ (^, 0), t^+00. (39) 



Thus, this solution describes the competition between two species when the species u even- 
tually dominate while the species v die. 

5. Conclusions 

In this paper, new definitions of Q-conditional symmetry for systems of PDEs are presented, 
which generalize the standard notation of non-classical (conditional) symmetry. It is shown 
that different types of Q-conditional symmetry generate a hierarchy of conditional symmetry 
operators. Since conditional symmetries can be applied for finding exact solutions of the relevant 
equations, which are not obtainable by the classical Lie method, we demonstrated this for 
constructing new exact solutions of the diffusive Lotka-Volterra system. 

Systems of DEs to find Q-conditional symmetries of two types for the nonlinear RD system 
([1]) are constructed. The case of conditional invariance under operator flT3|) with the coefficients 
(^ = 0, ril{t,x,U) and r]l{t,x,V) is analyzed in details. Using the recent paper [T3], we 
established that there are the nonlinear RD systems of five types, which possess Q-conditional 
symmetry (non-classical symmetry) operators satisfying definition 3. However, only one of 
them is simultaneously a Q-conditional symmetry operator of the first type. 

Definitions 1-3 can be straightforwardly extended on an arbitrary Tn-component system of 
PDEs presented in a 'canonical' form (the system has a simplest form and there are no any 
non-trivial differential consequences). However, new difficulties may arise because differential 
consequences of additional conditions 

Q{ui)=0,...,Q{um) = 0, (40) 

generated by the conditional symmetry operator ([6]). In fact, there are examples of single 
hyperbolic PDEs when the set of Q- conditional symmetry operators can be extended if one takes 
into account such differential consequences. For example, this occurs in the case of nonlinear 
hyperbolic equation Uu = uu^x- The reader may easily check that there is a much wider set of 
Q-conditional symmetry operators in [15] (Chapter 5) than it was found in [27] (Supplement 7) 
neglecting differential consequences. It is the reason why there is the definition of Q-conditional 
invariance (non-classical invariance) for a fc-th order single PDE, which requires to take into 
account all differential consequences of the additional condition Q{u) = up to the order k (see, 
e.g.,[2H|)- From this point of view, definitions 1-3 can be formally extended by using differential 
consequences of (1401) up to the order k = m.ax{ki, i = 1, . . . ,m}. However, we believe that 
such extensions produce much more cumbersome formulae but don't lead to any new operators 
for evolution system (jlj) because all the above mentioned differential consequences contain the 
mixed and/or high-order time derivatives, UtxiUttiUuxi ■ ■ ■■, which are absent in system (jl]). 
For example, if one applies the 'amended' definition 1 with Jvi\ = {Si = 0, S2 = 0, Q{u) = 
0,dtQ{u) = 0,dxQ{u) = 0} (instead of fITBl) ) to find Q-conditional symmetry operators of the 
first type for the RD system (fT2|) then the system of DEs (TT71) is exactly obtained but nothing 
more. 
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